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1. INTRODUCTION 

Meteorologists are inter- 
ested in modeling the vertical flow of 
heat and moisture through the soil in 
order to better simulate the vertical 
and temporal variations of the atmos- 
pheric boundary layer. In the present 
study the one-dimensional PBL model of 
Diertele (1979) is modified by the 
addition of transport equations to be 
solved by a finite difference technique 
to predict soil moisture. The model of 
Diertele is a refinement of the two- 
dimensional model of Bornstein (1975) 
and Bornstein and Roboc)c (1976) in which 
a one-dimensional time-dependent model 
of the planetary boundary layer simu- 
lates the vertical distribution of wind, 
temperature, and water vapor. 

2. FORMULATION 

Time dependent solutions for 
the vertical distribution of wind, 
temoerature, and water vapor in the 
planetary boundary layer (PBL) can be 
obtained from equations used by 
Bornstein (1975) in his URBMB7 (for ur- 
ban meteorology) model. This system of 
equations describes a two-dimensional 
(x-z plane) boundary layer model, with 
a lower, analytical "constant flux" 
layer of 25 meters, and an upper, numer- 
ical transition layer of 1900 meters. 

A one-dimensional version of 
the model valid over homogeneous terrain 
was developed by Diertele (1979) to in- 
clude the effects of: 1) heating due to 

divergence of the long and short wave 
radiative fluxes, and 2) an energy con- 
servation equations for the air-ground 
interface, including radiative, sensible, 
and latent heat fluxes. Use of this 
equation requires inclusion in the 
model of a subsurface layer for predic- 
tion of soil sensible heat flux. In the 
current formulation, the differential 


equations of soil heat and moisture 
fluxes are solved using m explicit 
fin'te difference technique to obtain 
soil temperature and moisture respec- 
tively. Instead of using a parametri- 
zation for solving surface latent heat 
flux, it is obtained more accurately 
from the calculations of surface mois- 
ture and temperature. 

a. Tr ansition Layer 

’ The equations governing at- 
mospheric variables in a one-dimensional, 
hydrostatic, Boussinesq transition layer 
were given by Bornstein (1975). With 
the additional assumptions discussed 
above, the 'liquations become 
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where all symbols are defined in Appendix. 

The verticity and stream function 
approach, used in the original model, is 
retained in the present one-dimensional 
study for convenience. The eddy diffu- 
sivities Ku, K||, and Ky (set equal to Kh) 
are specified using the interpolation 
formula of O'Brien (1970). Details of 
the formulation of radiative flux diver- 
gence which enters into the thermal 
energy equation (2-3) for PBL, can be 
found in Santhanam (1980) as can the 


datailB of tha •quatlona for th« eonatant 
flux aurfaea boundary layar. Only a 
cloudlaaa and pollutant fraa atmoaphara 
ia eonaidarad in tha praaant atudy, which 
tharafora raatrieta tha fanarality of tha 
modal. 


b. 


Soil Laver 


Tha flow of moiatura through 
aoil can ba tnrittan aa 

W*- ‘ ‘ 


whara tha capillary action tom (firat 
tarn on tha right hand aida of tha aqua- 
tion) ia yanarally largar than tha con- 
voction tern (niddlo tarn) and tha 
gravity tern Uaat tern) . Tha flow of 
heat through tha aoil can bo written aa 
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where the firat tem on the right hand 
aide of the equation repreaents aenaible 
heat flux and the aecond tem representa 
latent heat flux. All values of the 
moisture dependent diffusivities 
and D*) » as well as the hydraulic conduc- 
tivity (Kvj) and the modified themal 
conductivity (X') were obtained from 
Sasamori (1970) . One difference between 
current fomulation and that of Sas 2 unori 
(1970) is the present use of moisture de- 
pendent values of soil heat capacity and 
density. 


c. Grid and Boundary Conditions 

The one-dimensional version 
of the interlaced grid of Fromm (1964) 
computes vorticity and stream function 
values at grid locations vertically dis- 
placed by one-half of a grid interval 
from locations of the velocity components. 
Studies by Clark (1970) and Taylor and 
Delage (1971) discuss the superiority of 
variable grid spacing in achieving high 
resolution near the surface, and the grid 
netwcrk used in the present model pos- 
sesses such resolution. A complete list 
of boundary conditions used ^n the pres- 
ent simulation can be found in Santhanam 
(1980) . 


In order to obtain the fol- 
lowing equation for the time-varying 
surface temperature, a balance is assumed 
at the surface between the computed net 
radiant flux and convection in the atmos- 
phere, conduction into the soil, and 
evaporation: 

Or + - Obt - * Pm*^p“*®* 




In order to obtain the follow- 
ing similar equation for surface moisture, 
a balance ia aaaumad between the etmoa- 
pherie surface water vapor flux and soil 
surface moisture flux: 
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Equations (9) and (10) are partial dif- 
ferential equations with ttfo unknonma, 
i.e., moisture and temperature. Solu- 
tions of theae two equations are obtained 
by using a Newton-Ralphson double itera- 
tive technique (Carnahan at al., 1969). 

The specific humidity at 
ground can be obtained from 

qg - qa * (11) 

where qg is obtained from an expression 
given by Murray (1967) and hg ia obtained 
from equations from Philip (1957) and 
Nappo U975) . 

3 . SOLUTION 

At each time step, the 
following procedure is performed: 
determination of soil layer profiles; 
construction of constant flux profiles; 
determination of the new values of fi, q, 
n, and U at the surface and at the lowest 
finite determination of the new q, 
f, u, v,-and U pi.'^files using newest 
available values rZ all parameters. De- 
tails can be found in Santhanam (1980) . 

4 . RESULTS 

This simulation attempts to 
reproduce the temperature and moisture 
fields in the atmosi^eric and soil boun- 
dary layers. Initial conditions for 
these simulations were taken as those 
existing at the beginning of the 5th 
general observational period at O'Neill. 
Nebraska, 1230 LST on August 24, 1953 
(Lettau and Davidson, 1957). Numerical 
integration was carried out for a 48 
hour period. 

Computed soil moisture values 
near the surface (approximately 0.25 cm 
below surface) are compared with those 
of Sasamori (1970) in Fig. 1. Soil 
moisture displays a minimum value in the 
afternoon as the surface dries because 
evaporation is greater than the moisture 
flux from below. A maximum value is 
found after midnight. Wetness values 
decrease from 5.4 % at 1200 LST on the 
first day to 3.1% at 1200 LST on the 
next day, indicating a continuous deple- 
tion of soil moisture. The depletion 
arises as the moisture flux from below 
the surface is less than surface evapo- 
ration and as nighttime surface recon- 
densation is also less than daytime 
evaporation. In general, the present 
results are similar to those of Sasamori 
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raaults arc aimilar to thoaa of Saaaiaori 
(1970) . 

Computod aurfaca relativa 
humidity valuaa (Fig. 2) eompara wall 
with thoaa of Nappo (1975), who obtainad 
hia valuas by dividing obaarvad mixing 
ratio valuaa at 10 cm into tha atmoa- 
phara (which ha aaaumad aqual to that at 
tha aurfaca) by tha aaturation mixing 
ration. Howavar, hia apacific humidity 
valuaa tdiich wara obtainad by multiply- 
ing ralative humidity by aaturation 
apacific humidity (which waa not axplic- 
itly praaantad in hia original papar) do 
not natch tha praaant predictionn (Fig. 
3) in amplituda and phaaa. Howavar, tha 
praaant calculationa do agree aomewhat 
better with thoae of Saaamori (1970). 

Eatimated Bowen ration (of 
aennible flux to latent heat flux) from 
varioua modela are shown in Fig. 4. 

Thoae of Suomi (conputed from obaerva- 
tiona of net radiation and soil heat 
flux) ahow an unrealiatic jump at 1500 
LST due to a too small latent heat flux. 
In addition, his values beyond 18 hours 
of simulated time, as well aa those of 
Lettau (Lettau and Davidson, 1957) , 
which are obtained in a similar manner 
to those of Suomi, do not agree with the 
present results and those of Sasamori 
(1970) . These later results aa well as 
the present ones, are from complete PBL 
formulations, while those of Suomi and 
Lettau are associated with simple para- 
metrizations. 

5. CONCLUSION 

The two dimensional URBMET 
urban planetary boundary layer model has 
been simplified to one dimension and 
modified to perform a detailed analysis 
of the flow of soil moisture into the 
atmosphere. Model calculations of tem- 
perature and moisture in the soil and 
atmospheric boundary layers were vali- 
dated against existing observations at 
O'Neill, Nebraska. 

Results have indicated that 
simple models with oarametrizations of 
the surface latent heat flux do well in 
predicting surface relative humidity but 
do not do well for specific humidity, 
because they do well in predicting sur- 
face temperature. Calculations show 
that models using empirical techniques 
to evaluate the surface energy balance 
do not predict Bowen ratios as well as 
those incorporating a soil layer. 
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Variabltt* and Constants - Croak Alphsbat 


APPENDIX 
List of Symbols 

Varisblss and Constants • Roman Alphabot 



soil hast capacity 

specific hast of dry air 

total soil moistura dif> 
fusivity 

vapor moistura diffusiv- 
ity 

total soil thermal dif- 
fusivity 

Coriolis parameter 

surface relative humidity 

eddy exchange coefficienti 
for heatf mcmicnttBii, and 
water vapor 

hydraulic condictivity 


latent heat of vaporiza- 
tion 

<]« g > <1 f specific humidity • sur- 
^ * face q, saturated q, 

and friction q 

Q , 0 vertical flu^ of latent 

^ " and sensible heats 

0 ,0 . upward and dovmward di- 

rected long-wave radia- 
tion 

0 , Q net all-wave and solar 

** " radiation at the 

suirface 


c 


time 


T 




surface and soil tempera- 
ture 


u, 


horizontal component of 
the wind and geostror-i'? 
wind in x-direction 


u* friction velocity 

V horizontal component of 

the wind in the y- 
direction 

X, y, z horizontal coordinate iu 

direction of geos tropic 
wind, perpendicular to 
geostropic wind and 
vertical coordinate 


S • 8* 
X' 


o 

T 


wave length integrated 
surface emisslvlty 

modified y-cmaponent of 
vorticity 

soil wetness 

potential temperature end 
friction • 

modified soil, thermal 
conductivity 

constant space averaged 
atmospheric density 

density of water 

Stefan-Boltzman constant 

stream function 



Predicted soil moisture (by volume) at 
0.25 cm below surface, where 12, the 
initial time, corresponds to 1200 LST. 








